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The proton light output function in electron-equivalent energy of various scintillators based on 
linear alkylbenzene (LAB) has been measured in the energy range from 1 MeV to 17.15 MeV for 
the first time. The measurement was performed at the Physikalisch-Technische Bundesanstalt 
(PTB) using a neutron beam with continuous energy distribution. The proton light output data is 
extracted from proton recoil spectra originating from neutron-proton scattering in the scintillator. 
The functional behavior of the proton light output is described succesfully by Birks' law with a 
Birks constant kB between (0.0094 ± 0.0002) cmMeV"^ and (0.0098 ± 0.0003) cmMeV"^ for the 
different LAB solutions. The constant C, parameterizing the quadratic term in the generalized 
Birks law, is consistent with zero for all investigated scintillators with an upper limit (95% CL) of 
^ ' about 10~^ cm'^ MeV~^. The resulting quenching factors are especially important for future planned 

supernova neutrino detection based on the elastic scattering of neutrinos on protons. The impact 
of proton quenching on the supernova event yield from neutrino-proton scattering is discussed. 
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I. INTRODUCTION 



' In the last decades enormous progress has occured in the field of neutrino physics. A new field of astrophysics 
Qh, has been opened with the observation of solar and supernova neutrinos [H-Q • First indications of ultra-high energy 
Q ' neutrino detection exist Q . From all sources the appearance of supernova neutrinos is the least predictable but a lot 

\ can be learned if such an event occurs in the Milky Way. In the final stages of a core-collapse of massive stars nearly all 

■ the binding energy of the star is expected to be released in the form of neutrinos and antineutrinos of all flavors. The 
' expected spectrum consists basically of two components: The deleptonization burst consisting of released within 

a few ms and the emission of all kinds of neutrinos from the Helmholtz-Kelvin cooling phase of the protoneutron 
star lasting several seconds. For recent reviews see [^-[ll- The intense neutrino pulse can be observed on the Earth 
J> , by neutrino detectors with a sufficiently low energy threshold of about 1 MeV. For the SN 1987A this was achieved 
' by Kamiokande II, 1MB and the Baksan Scintillator Telescope 041] ■ The signal is expected to be dominated by Ve 
\ because of the large cross section for the respective reaction. However, it is of great interest to measure the fiux of 

■ Vr, Vfi and collectively called v^, and as well. The detailed shape of the overall spectrum depends strongly on 
the model simulation. Recent progress has been made in using 3-dimensional simulations taking into account large 
scale convection. Furthermore magnetic fields, general relativity and other effects were included resulting in different 

' spectral details. As the average neutrino energy is below about 30 MeV, charged current interactions are only possible 
I for electron type neutrinos, leaving only the far more challenging neutral current reactions as possible experimental 
detection channel for supernova Vx- In the references [^, neutrino-proton elastic scattering, v + p ^ v -\- p, \n 
liquid scintillator detectors is proposed for the measurement of Vx since it is the only neutral current channel providing 
spectral information. Though the total cross section of this process is about a factor of four smaller than the cross 
^ ^ section of inverse-beta decay, the yield above a realistic threshold of about 200 keV is of the same order since this 
d • reaction is possible for all six neutrino types foj. 

This detection channel plays an essential role for large scale deep underground liquid scintillator detectors like 
KamLAND, Borexino and SNO+, in the pending measurement of the total energy and temperature of v^, Vr, 
and i>r- It should be noted, however, that only a fraction of the deposited kinetic energy of highly ionizing particles 
like protons is visible. Hence to properly reconstruct the true kinetic proton energy Ep from the visible energy E^^^, 
the light response of the liquid scintillator for these particles has to be known. The ionization quenching of the light 
yield is well-described by Birks' law 11] and can be quantified in terms of the parameter product kB known as 
Birks' constant and a second parameter C. While the proton light response function of the KamLAND scintillator, 
which consists mainly of pseudocumene (80%) and normal-paraffine (20%), is already measured up to 10 MeV [lH, 
the response to protons of scintillators based on linear alkylbenzene (LAB) has not been measured before. LAB is 
employed for novel liquid scintillators which are used or planned as the neutrino target in recently commissioned. 
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upcoming and potential future neutrino experiments like Daya Bay, RENO, SNO+ and LENA. Also further potential 
experiments like LENS, HANOHANO and Daya Bay II are considering LAB as possible solvent for their scintillator. 

This article presents the measurement of the light response function L{E) for protons relative to the electron 
response function Le{E), carried out at the accelerator facility PIAF of the Physikalisch-Technische Bundesanstalt 
(PTB) for different LAB based scintillators. Recoil protons are produced in the liquid scintillator by neutron-proton 
scattering events using a neutron beam with continuous energy distribution. The maximum energy of the recoil 
protons is determined from the energy of the incident neutrons, which is measured using the time-of-flight (TOE) 
method. Hence, the proton light output function can be determined from the position of the recoil proton edge in the 
pulse-height spectra produced by mono-energetic neutrons. The presented method was already used successfully at 
the PTB for the characterization of NE21^1, BC5O10 and BC501A scintillation detectors [H-Ell. 

II. MEASUREMENTS 

In the present experiment, the light output function of binary and ternary LAB based scintillation systems is 
determined. In these solutions, 2,5-diphenyloxazole (PPO) is employed as primary fluor and p-bis-(o-methylstyryl)- 
benzene (bis-MSB) as secondary fluor, which acts as wavelength shifter. Two concentrations of PPO in the LAB 
solvent, 2g/l and 3g/l, are studied, each with and without bis-MSB, for which a concentration of 15mg/l is chosen. 

After dissolving the solute in LAB, the solutions were purified from oxygen by bubbling with argon for 30min 
before filling the samples, without further contact to air, in a scintillation detector consisting of a cylindrical dural 
cell with one port covered by a window made of ground and polished fused silica for higher UV transparency. The 
diameter and inner height of the cell is 50.8 mm. The cell is coupled to a XP2020Cjl photomultiplier tube (PMT), 
also equipped with a fused silica window, by a conical UV transparent poly(methyl methacrylate) (PMMA) light 
guido The inner walls of the cell and the light guide are coated with BC-622i^ reflective paint to increase the 
light collection efficiency. The detector is equipped with a LED gain stabilization system which regulates the high 
voltages such that the maximum drift of the gain is smaller than 1%. Standard NIM modules are used to obtain a 
pulse-height (PH) signal from the ninth dynode out of twelve, as well as a pulse-shape (PS) signal and a time-of- 
flight (TOE) signal from the anode of the photomultiplier. The ninth dynode is chosen for obtaining the PH signal 
to avoid a non-linear gain of the PMT. The PH, PS and TOE signals are registered by a multi-parameter data 
acquisition system and stored event -by-event. The PS discrimination module, as described in detail in JJj, uses 
the zero crossing method to derive a timing marker from the trailing edge of the anode pulse. The time difference 
of this PS timing marker to the TOE timing marker derived from the fast leading edge of the anode pulse, using 
a constant fraction discriminator (CED), is measured with a time-to-amplitude converter (TAG). It is related to 
the decay time of the scintillation light and therefore to the ionization density of the charged particle producing the 
event. The non-linearity of the TAG converter is measured using a time calibration module which produces a set 
of pulses at multiple intervals of At = 40 ns. The variations of the peak positions from the nominal values quantify 
the non-linearity and are determined to be smaller than ± 0.23 ns. Eig. [T] shows the registered events sorted into 
a PS versus PH matrix independent of the energy of the incident neutron. The structures due to neutron-induced 
events ('n'), like recoil protons or other light charged particles, are partially separated from those caused by Gompton 
electrons ('7'). The electron events are due to ambient 7-radiation and photons resulting from the inelastic neutron 
scattering on carbon nuclei in the scintillator and on alumnium nuclei of the walls. A suppression of events with a PS 
signal below the indicated line strongly reduces the influence of photon-induced events with respect to the PH spectra 
of neutron-induced events. However, for LAB based scintillators a fraction of misidentified photon-induced events 
remains in the overlap region, as indicated by Eig. [T] This is mainly effected by the small number of photoelectrons 
produced per scintillation event which is primarily due to the wavelength mismatch between the emission spectrum 
and the sensitivity of the photocathode. 

The energy calibration of the PH spectra and determination of the light response to electrons for each scintillator 
is obtained from the position of the Gompton edges produced by a set of standard 7-ray sources (^^^Gs, ^^Na, ^°^Bi). 
Ambient background is measured separately and subtracted from each source measurement. The analysis of the 
calibration data is described in detail in Sec. IIIII 

To obtain a proton light output function over an energy range as large as possible in a single experiment, a neutron 
beam with a continuous spectral distribution is used, which is produced at the PTB TOE facility by bombarding 
a stopping-length beryllium target with a 19MeV proton beam from the isochronous GV28 cyclotron [3]. At an 
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emission angle of 0°, the maximum energy of the neutrons is 17.15 MeV. The center of the ceU was positioned at a 
distance of 12.11m from the Be target under 0° with respect to the proton beam. The distance between target and 
detector cell is measured with an accuracy of ± 2 mm. To consider the target thickness of 3 mm and include the 
uncertainty of the effective center of the cell caused by the flucnce attenuation, ± 6 mm are conservatively assumed. 
With a beam current of about 50 uA, the neutron rate at the detector front face was about 3 x 10^ s~^. The repetition 
frequency of the proton beam was set to 481.3kHz. The trigger rate of the CFD module was about 3 x 10'^ s^^ for a 
threshold set close to the electronic noise. The probability for multiple neutron interactions within one beam pulse 
is thus smaller than 1% and does not cause a deterioration of the shape of the pulse-height spectra. In the given 
configuration the neutron energy threshold for time-frame overlap is 0.15 MeV, which corresponds to an electron- 
equivalent energy for the given scintillators of less than 20 keV. Hence, with a hardware threshold of about 200 keV, 
time-frame overlap is avoided. 

Proton bunches that are not entirely deflected by the internal beam pulse selector system cause additional pulses, 
so-called satellite pulses, on the target. The spectral distributions of outgoing satellite neutrons are assumed to be 
the same as the ones from the main proton bunches and are visible in the PH vs. TOF matrix as faint distributions 
which are shifted in TOF relative to the distribution resulting from the main beam pulse. Hence, satellites cause 
spurious events above the recoil proton edge in the pulse-height spectra which could affect the determination of the 
position of the recoil proton edge. A subtraction of the shifted TOF spectra, though, showed no significant influence 
on the position determination of the proton recoil edge. 

To improve the dynamic range of the pulse-height measurement for neutron energies below 5 MeV, each measure- 
ment is repeated with an increased amplification. This amplification mode is referred to as high gain (HG) mode, the 
other one as low gain (LG) mode. For the HG measurement, n/7-separation by pulse-shape discrimination (PSD) 
is not possible. However, at these energies, inelastic scattering on carbon does not contribute to the neutron induced 
response because of the energy of the first excited state in ^^C which is at 4.439 MeV. The contribution from prompt 
7-radiation from the target and from ambient 7-radiation is subtracted experimentally. 



III. ANALYSIS 



1. Calibration 



The amount of detected scintillation light and the pulse-height resolution strongly depend on the specific composi- 
tion of the scintillator. Therefore the pulse-height calibration of the scintillation detector is repeated for each filling 
and amplification mode using three 7-ray sources with a total of six 7-lines, namely ^'^^Cs [E^ — 0.662 MeV), ^^Na 
{E-y = 0.511 MeV and E-^ = 1.275 MeV) and ^o^Bi [E^ = 0.570 MeV, E^ = 1.064 MeV and E^ = 1.770 MeV). 

The position of each Compton edge, with energy 



0.bll + 2E^' ^ ' 

where all energies are given in MeV, is defined by comparing the measured pulse-height spectrum with a distri- 
bution simulated using the Monte Carlo photon transport code GRESP [l^. To adjust the theoretical shape to the 
measurement, the simulated spectrum is folded using a Gaussian resolution function with a pulse-height dependent 
width and fitted to the measured pulse-height spectra (Fig. [5]). 

In this way the measured pulse-height signal PH is calibrated in units of the maximum kinetic energy E,, of the 
Compton electrons. The calibrated pulse-height (light output) L is given by 

L = n- PH, (2) 
where n denotes the calibration factor. The relative resolution AL at pulse-height L is parameterized by [lO] 



The fit parameters a, /3, 7 can be attributed to the individual contributions from spatial dependence of the light 
collection efficiency, statistical variation of the number of photoelectrons and electronic noise, respectively. The values 
of the resolution parameters for each scintillator filling are obtained by fitting Eq. ([3]) to the locally determined full 
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width at half maximum AL at each Compton edge position. Fig. [3] shows the resulting resolution functions for every 
sample as a function of the light output. The lower resolution for the scintillators without wavelength shifter is mainly 
attributed to detector properties. The reflectivity of the BC-622A coating inside the cell is reduced below 420 nm 
and the XP2020Q PMT has its peak sensitivity at 420 nm, while the emission peak of PPO lies around 360 nm. Thus 
the spatial dependency and statistical variation are increased. However, the properties are more fortunate for the 
scintillators with bis-MSB which has a maximum emission around 420 nm. For this reason, besides the two binary 
scintillators, the two systems with additional wavelength shifter are investigated, though no difference in the ionization 
quenching properties is to be expected. 

According to the light output function L{E) is related to the stopping power dE/dx of a charged particle of 
kinetic energy E in the traversed scintillator via 

L{E) ^ S - [ de 
Jo 

if the particle is stopped in the scintillator. This relation is known as Birks' formula with the product kB being 
the Birks constant. The light output L is the total light emitted when a charged particle loses all of its energy E 
within the scintillator. For electrons with energies > 125 keV, the stopping power becomes very small jllj and Eq. (|4]) 
reduces to the linear expression 

Le{E) ^S-iE^ Eo). (5) 

The subscript 'e' is added here to underline that this description of the light output only holds for electrons. The 
parameter S denotes the chosen scaling between electron energy and light output and is set to 5 = 1, which yields 
the light output in electron-equivalent energy, if the calibration procedure described above is used. The constant Eq 
accounts for the fact that the response of the scintillator to electrons is only approximately linear. The non-linearity 
at small energies leads to an energy offset, which is assumed to be Eq — 5keV for NE213 and BC501. Recent 
measurements with LAB and EJ-301 scintillators [2l|, though, revealed a non-linear behavior already starting at 
400 keV, resulting in an offset of about 50 keV. Within the present work, the potential non-linearity of the electron 
light output can not be discriminated from the effect of the remaining electronic offset of the amplifiers which was 
corrected for during the analysis of the pulse-height calibration, i.e. Eq is set to zero. 

A comparison of Eq. ^ and Eq. ([5]) gives the relation between PH and electron energy 

PH = n-^ ■ {E ~ Eo) 

= n-^ ■ E, (6) 

where the second equality only holds for the case described here. The calibration factor is now derived from a 
linear fit to the positions of all Compton edges in the PH spectrum as performed in Fig. 2] and is independent of the 
total offset of the pulse-height scale. The values obtained in the present measurements are summarized in Tab. HI 

2. Proton light output function 

For the determination of the proton light output function L{E), basically the same technique is used as for the 
analysis of the photon-induced events in the calibration measurements described in Sec. IHI II The events measured 
with the white neutron beam are first sorted into a TOE vs. PH matrix. As mentioned above, PSD is used to select 
only events produced by light charged particles for the spectra obtained with the LG setting of the amplifier, while in 
the HG mode n/7-separation by PSD is not feasible. The TOE range covers the full neutron energy range from about 
IMeV to 17.15MeV. The TOE scale is established using the arrival time difference between prompt gammas and 
neutrons of maximum energy and a time determined using an electronic time calibration module. The prompt gamma 
line is broadened due to the finite time resolution of the accelerator and the uncertainty of the centroid position of 
the resulting peak is calculated to be ± 0.1ns. Calibration gives a TOE channel width of 0.7764ns and the relative 
uncertainty of the time calibration is 0.05%. Pulse-height spectra induced by quasi mono-energetic neutrons are then 
extracted by collecting all events contained in small windows around the TOE value corresponding to an individual 
neutron energy. The size of the TOE window is always made smaller than the pulse-height resolution AL for protons 
of the respective energy. A second window of the same size is placed on the random background at TOE values 
smaller than those of neutrons with the maximum energy and the extracted background spectrum is subtracted from 
the signal spectrum. Pulse-height spectra for different incoming neutron energies are shown in Fig. [5] The energy of 
the recoil protons, producing the edge in the pulse-height spectra, equals that of the incident neutrons. 
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Before extracting the position of the recoil edge, the dependence of the TOF measurement on the pulse-height is 
investigated. The so-called time-walk is caused by imperfections in the constant fraction timing technique employed 
in the CFD module which show up in particular at small pulse-height. It leads to a distortion of the PH vs. TOF 
matrices and finally in the extracted proton recoil spectra at low pulse-heights. A correction of the PH vs. TOF 
matrices of LAB with 2 g/1 PPO and 15 mg/1 bis-MSB revealed a deviation of about only one channel in the TOF 
spectrum. For comparison, the smallest window width, which is chosen at high neutron energies, has a width of five 
channels. A spectral distortion in the region of small pulse-heights, however, has the biggest infiuence on the definition 
of the proton recoil edge at low neutron energies, where the TOF window width has to be increased because of a 
decreased resolution, as described in Sec. IIII 21 The maximum width of the TOF windows amounts to 25 channels, 
which is much larger than the observed deviation in the measured TOF. Hence, time-walk has a minor effect. As 
a correction of the PH vs. TOF matrices, to account for time-walk, can bring in new artefacts due to the binned 
structure of the data, no correction is performed within this analysis. 

To calculate the detector response to monoenergetic neutrons, the Monte Carlo code nrespT [2^ is used. nrespT 
models all reactions of neutrons in the scintillator, in the detector housing and in the Lucite light guide, including the 
production of secondary charged particles, but it does not simulate the interaction of deexcitation photons (e.g. from 
the first excited state in ^^C at 4.439 MeV) resulting from inelastic neutron scattering. Instead, it is assumed that such 
events are suppressed by PSD. If the incoming neutrons have enough energy to excite ^^C, signals from deexcitation 
photons with pulse-heights up to about 180 channels add up to the proton recoil spectra. As this background is 
not fully suppressed by PSD, it leads to more counts above threshold in the measured pulse-height spectra than in 
the simulated ones (see Fig. for neutron energies higher than ^4.5 MeV. However, this mainly affects the proton 
recoil edge determination, if the recoil edge is in the same region as the photon induced background. This is only 
the case for a maximum proton energy (i.e. incoming neutron energy) up to ^^TMeV. For the fit of calculated to the 
measured spectrum, a pulse-height interval is chosen, which starts below the proton recoil edge and ends above. For 
neutron energies above about 7 MeV, this interval is chosen such, that photon induced background is excluded. The 
pulse-height given in channels is finally translated into light output in electron-equivalent energy L using the factor 
n^^ obtained in the calibration with photon sources. The resulting relative proton light output resolution functions 
are plotted in Fig. [6l 

The simulated pulse-height spectra are adapted to the experimental ones in an iterative process. As an initial 
step, an existing reference light output function for NE213 detectors is used to calculate pulse-height spectra using 
nrespT. The simulated spectra are compressed or stretched and adjusted in height to fit the experimental spectra 
in the pulse-height region around the recoil edge and the respective compression or stretching factors are applied to 
correct the light output function. The corrected function is inserted as input for the next iteration step. The behavior 
of the applied light output function follows a straight line for > 8 MeV, where the parameters of the line are gained 
from a fit to the extracted data points. After about four steps, the procedure converges. Fig. [7] compares the initial 
light response function with the final one for LAB with 2 g/1 PPO and 15 mg/1 bis-MSB. The data points extracted 
in this way are presented in Fig. [S] together with the electron data points gained from calibration. The results for the 
remaining scintillators are shown in Fig.lHl- fTTl The contributions to the total uncertainty in L and E are summarized 
in Tab. IB 

The most popular analytical description for the energy dependency of the light output is given by Birks' formula 
Eq. (|4]), assuming ionization quenching to be the reason for the reduced light output. This expression, though widely 
used for all kinds of scintillators, was originally developed for anthracene crystals. However, the underlying physical 
processes in liquid scintillators are different and this description should be considered as being semi-empirical. For 
large values of dE / dx, like for incident ions, Eq. Q is often encountered in its generalized form 



L{E) = S- 
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1 + kB 



dE 
dx 



C 



dE 
dx 



(7) 



including a quadratic correction term parameterized by C, as proposed in [23j . Within the present work, electron 
and proton light output functions of a particular scintillator are measured under identical conditions so that the 
scaling factor S for proton data is the same as the one for electron data, i.e. = 1. The Birks constant kB and the 
parameter C for LAB based scintillators are obtained by means of a fit of the theoretically expected light output 
^thoo Q ^Yie experimentally determined values L°^p presented in Fig.|8]-[TT] Extracted data points, instead 

of the final input function in nrespT, are used for this fit to be independent from the mentioned straight line fit for 
E > 8 MeV. The stopping powers dE/dx for protons are calculated using the code SRIM [23|. The calculation is 
carried out following the pull approach [251] . where 
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Within this calculation K = 8 systematic uncertainties, listed in Tab. HH are included as nuisance parameters. 
describes the fractional change of the n-th value of £*heo j-f |-j^g fc-th source of systematics is varied by 1 dk and £,k is 
a standard normal deviate. The normalization condition for the ^fc's is realized through quadratic penalties, summed 
over the K sources of systematics. The statistical uncertainty, denoted by w, is also indicated in Tab. [TTl For each 
combination of kB and C, Eq. ^ is minimized with respect to all ^fc's using minuit and the minimal value, Xpuii' 
stored in a {kB, C) map. The best fit values are finally extracted at the minimal XpuW ^^1^ parameter space. The 

uncertainties are extracted from the Xpuii projections on the kB resp. C axis. The best fit value of kB has a two-sided 
limit and the kB values at Xpun — Xpuii min = ^ S^'^^ ^'^ uncertainty. C has a one-sided limit and the C value at 
Xpui\ ~ Xpuii min ~ 1-645 rcfcrs to an upper limit at 95% confidence level. Besides the experimental uncertainties, a 
further uncertainty has to be considered, namely the uncertainty of the stopping power calculation entering Eq. ([7]). 
With the SRIM code, a determination of the stopping power for elemental materials is possible with an accuracy of a few 
percent [l^, thus the described fitting procedure is repeated for dE/dx + adE/dx and dE/dx — OdE/dx- Assuming a shift 
of 2% of the stopping power already leads to a 2% difference in the resulting quenching parameters. This uncertainty 
is added quadratically to the experimental one and is the dominant contribution to the resulting total uncertainty. 
The fit results are summarized in Tab. IIIII C is consistent with zero for all investigated scintillator samples with an 
upper limit (95% CL) of about 10"'^ cm^ MeV"^ and the values for kB range between (0.0094 ± 0.0002) cm MeV"^ 
and (0.0098 ± 0.0003) cm MeV"!. 

According to current knowledge, all ionization quenching processes are primary processes in the scintillator, i.e. 
processes that transfer ionization energy to excitation energy of the solvent. Quenching processes compete with the 
excitation of the solvent molecules into 7r-electron singlet states. This is the primary process which is essential for the 
final scintillation emission in unitary systems as well as in binary and ternary systems. Consequently, for scintillators 
with the same solvent, the magnitude of ionization quenching should be the same independent of possible solutes. 
Within the presented measurement no significant deviation from these considerations could be observed. 

For practical use, the proton quenching factor Qp, the ratio of proton and electron light output functions 



QAB) = Ifi. (9) 

for each examined LAB scintillator is presented in Fig. [121 Above a kinetic energy of 8 MeV, the proton light output 
is about half the electron light output. Below that energy the quenching factor strongly decreases, i.e. the proton 
light output is further reduced down to 20% of the electron light output for i? w 1 MeV. 



IV. APPLICATION OF THE RESULTS 



The detection of supernova neutrinos of all fiavors using neutrino-proton elastic scattering is discussed in detail in 
[lo| . Within that reference, different experiments capable of the respective detection are compared, namely Borexino, 
KamLAND, SNO+ and LENA. For the LAB scintillator detector SNO+, a Birks' constant of kB = 0.0073 cm MeV^^ 
is used. However at that time no measurement of proton quenching in LAB was published. The measurements 
presented here for LAB with 2g/l PPO (as will be used in SNO+) reveals with kB = 0.0098 cm MeV~^ a stronger 
proton quenching. To show the impact of a stronger quenching on the event yield induced by a supernova, the same 
detector (i.e. SNO+) and assumptions for the supernova are used as in In the same reference, the prospect of the 
larger detector LENA is discussed, assuming kB = 0.01 cmMeV"^ and 3.3 x 10^'^ free protons in the liquid scintillator. 
As LENA favors LAB as solvent with 3 g/1 PPO and 20 mg/1 bis-MSB "^T^ , this detector is also discussed in the present 
article with the obtained result of kB = 0.0098 cm MeV~i for LAB with 3 g/1 PPO and 15 mg/1 bis-MSB. For both 
detectors, C is considered to be zero in agreement with the results summarized in Tab. IIIII 

The supernova is assumed to occur at a distance d = lOkpc, releasing over the duration of the burst a total energy 
of e = 3 X 10^'^ erg which is equipartitioned among all flavors. The energy is distributed according to 
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— = (10) 

dE^ ^ dE^ ^ ^ 

a 

_ 2.35 X 10" y-^Ea E^ / AE^ \ 
~ cm^MeV ^ '^Je^J'''^ ' 

where dFa/dEi, is the neutrino fluence per flavor. The energies are given in MeV, d in lOkpc and Eq, in 10^^ erg. 
The mean energy {E^)a of the individual flavors differs and is exemplarily set to 12 MeV for v^, 15 MeV for and 
18 MeV for v^, foflowing the approach in [l3|. The respective fluence in v^, Ve and the sum of afl as well as the 
total fluence are shown in Fig. [T31 The fluences are at the order of 10^°cm~^ MeV~^ and the neutrino energy E^ 
ranges up to about 60 MeV. 

For these low neutrino energies, the differential neutrino-proton elastic scattering cross section [2^ can be simplified 

to 



da _ Gp 

Ie'p ^ 



TT 



4.83xlO--cmV^_^4gg^) (11) 



MeV V 

(considering only the lowest order in Ey/rup) and the cross section for neutrinos and anti-neutrinos can be regarded 
as identical 9]. For the last step, a proton mass of nip = 938 MeV is used and the neutral-current coupling constants 
are cy = 0-04 and ca = 0.635 [9]. Ey and the proton recoil energy Ep are given in MeV. 

Using Eq. (|10p and Eq. (jlip . the true recoil energy spectrum of recoil protons produced by supernova neutrinos is 
given by 

— -nH dE —— (12) 
dEp J dEu dEp 

where Np is the number of free protons in the scintillator target. The minimum neutrino energy needed to accelerate 
a proton to energy Ep is tv^min ~ \J mpEp/2 [gj. Fig. [T4l shows the expected proton recoil spectra for the various 
flavors in a detector with Np = 5.9 x 10'^^ (corresponding to the SNO+ target mass of about 0.8ktons [1^). Despite 
the small cross section, the intense neutrino fluence apparently causes a significant neutrino-proton elastic scattering 
yield is of the order of 10^ events per kton LAB. The proton energies ranges from OMeV to about 5 MeV. 

However, proton quenching markedly changes the shape of proton energy spectra and shifts the true proton recoil 
energies to lower observed energies. Even though there is no energy threshold to produce scintillation light, below 
0.2 MeV organic liquid scintillators are strongly affected by /3-decays of ^'^C, which is intrinsic to organic solvents. 
The high decay rate dominates this region and thus sets an effective threshold of 0.2 MeV (further backgrounds are 
not taken into account within this article as they are described in detail in '9^1 and found to be negligible for the 
presented supernova signal channel). Hence, the observed event yield finally depends on the quenching strength. The 
expected event yield above threshold will be discussed in the following for SNO+ and LENA. The respective detector 
specifications are given in Tab. II VI The main difference between both detectors is the amount of target protons, which 
is about 5.9 x lO^^ in SNO+ and 3.25 x lO^^ in LENA, given a proton density in LAB of 6.31 x lO"^"^ cm.'^ [sO]. 

To convert the true proton energy spectra in Fig. [Uinto spectra of visible energy Ep^^ in electron equivalent, the 
quenching factor Qp{Ep), as presented in Fig. [121 is applied to Eq. p2|) . The resulting spectra are folded with the 
energy resolution and scaled according to the respective number of free protons as given in Tab. II VI Fig.ll5land Fig. [TCI 
finally show the expected supernova neutrino signal spectra from v-p elastic scattering in SNO+ and LENA under the 
given assumptions for the supernova scenario and the detectors. The event yield above a threshold of 0.2 MeV, using 
kB = 0.0098 cm MeV"i and C equal to zero for LAB, is about 98 events in SNO+ and 5403 events in LENA. The 
yield scales linearly with the target size. For LAB scintillators it is in the range of ~123 events/kton, which makes 
this detection channel very promising for large scale detectors. Applying the listed energy resolution only changes 
both yields by less than 0.5% compared to perfect resolution and seems to have a negligible effect. The dependence 
on energy resolution thus is not further studied here, however needs to be investigated when a profound sensitivity 
study to supernova neutrino fluence parameters is performed. 

As mentioned earlier, the determination of kB is most sensitive to the calculation of the stopping power dE/dx 
and vice versa the calculation of Qp for a fixed Birks' constant. In reference ^Ifl] dE/dx tables from pstar [s^I are 
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used while for the presented analysis SRIM tables are employed. Though the stopping power data is similar, the 
expected yield for SNO+ is recalculated with kB — 0.0073 cm MeV~^ and C considered to be zero, using SRiM tables, 
to show the impact of a higher kB value without bias. Fig. [IT] compares the true proton energy sum spectrum and 
the observed sum spectra assuming kB — 0.0073 cm MeV~^ and kB — 0.0098 cm MeV~^ respectively. The higher 
quenching reduces the expected event yield by about 16% from ~116 event^ to ^98 events. This figure demonstrates 
how proton quenching pushes supernova v-p elastic scattering signal events to lower energies and that the yield 
above threshold strongly depends on the strength of the quenching. Therefor it is of great importance for liquid 
scintillator detectors, sensitive to supernova neutrinos, to properly measure the proton quenching parameters of the 
used scintillator to extract valuable information about supernova neutrino fluences of all flavors. 

V. SUMMARY AND OUTLOOK 

This article presents a measurement of the light output function L{E) of protons relative to the one of electrons for 
four different LAB based scintillators. A measurement of the relative proton light output reveals the magnitude of light 
output reduction due to ionization quenching. Fitting the functional description of L(E) proposed by J. B. Birks, pa- 
rameterized by kB and C, to the experimental data yielded values for kB ranging between (0.0094 ± 0.0002) cm MeV~^ 
and (0.0098 ± 0.0003) cm MeV-^ and an upper limit of about IQ-'^cm^ MeV'^ (95% CL) for the parameter C for 
the investigated scintillators. 

It should be noted, however, that the parameters within Birks' law, kB and C, alone are not conclusive for the 
comparison of different scintillators (i.e. scintillators with different solvents) as they are correlated with the stopping 
power of the scintillator material. Consequently, to decide between two different scintillators, always the full light 
output functions need to be compared. Different solutes in the same solvent, on the other hand, do not show a 
difference in the calculated stopping powers due to their small concentration. For the latter case, kB and C are 
thus sensitive to an eventual influence on the strength of ionization quenching by different solutes. The presented 
measurement does not reveal a significant disparity of the resulting parameters of all LAB based scintillators and is 
in accordance with the expectation that all ionization quenching processes are primary processes. 

A further remark is that the values for kB extracted from measurements of the non-linearity in the light output 
of low-energy electrons is not necessarily identical with the one of heavy ionizing particles like ions. Whereas it 
is assumed that the kB value is the same for different ions, like protons and alphas, once the conditions of the 
measurement are fixed. This hypothesis can be tested with the taken data by determining the a light output function 
for the investigated scintillators. ^^C(n,a)^Be and ^^C(n,n'3a) reactions, induced by neutrons from the beam, lead 
to characteristic structures in the pulse-height response of the scintillator which can be reproduced with the nresp7 
code. A proper knowledge of the scintillator response to a's is important for the development of a sophisticated 
background model of a liquid scintillator detector, including a-induced backgrounds. This analysis and discussion 
though is beyond the scope of this article and an outlook to future possibilities. 

Using the measured proton light output functions for LAB solutions, the potential of LAB based scintillator detectors 
to observe supernova neutrinos of all flavors via v-p elastic scattering is discussed and it is shown how a different 
Birks constant affects the event yield. Due to this impact on the observation, well-known quenching parameters are 
important for liquid scintillator detectors, capable of supernova neutrino detection. 
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Table I. Calibration factor n ^ of the investigated LAB based scintillators obtained in LG and HG mode. The given la 
uncertainties are the results of a least-squares fit to the data points shown in Fig. |4] 



LAB admixture 


Gain 








[MeV^^] 


2g/l PPO, 15mg/l bis-MSB 


LG 


324.5 ± 1.0 




HG 


343.4 ± 1.1 


2g/l PPO 


LG 


297.8 ± 1.4 




HG 


324.0 ± 2.1 


3g/l PPO, 15mg/l bis-MSB 


LG 


325.5 ± 0.9 




HG 


350.5 ± 0.7 


3g/l PPO 


LG 


294.7 ± 2.0 




HG 


325.1 ± 2.5 



Table II. Systematic and statistical uncertainties of the measurement of the relative proton light output L{E) and their la 
value. Time-walk of the CFD, satellite pulses, multiple neutron events and time-frame overlap have a minor effect and are 
neglected. The extracted proton recoil edge position carries two uncertainties. The statistical uncertainty is the uncertainty on 
the compression factor within the adaption of calculated to measured proton recoil spectrum. Furthermore, since the choice 
of the intervall for the fit shifts the recoil edge position by more than IcTstat., this shift is considered as additional systematic 
uncertainty. 



Systematic uncertainty 


1 CTsyst. 


prompt 7 peak centroid position 


±0.1 ns 


TAG non-linearity 


±0.23 ns 


time calibration 


±0.05% 


target - detector distance 


±6 mm 


gain stabilization 


±1% 


pulse-height offset 


±2ch 


recoil edge position 


±2% 


calibration factor 


(see Tab. HJ 


Statistical uncertainty 


1 0"stat. 


recoil edge position 


<0.02% 
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Table III. Quenching parameters kB and C, following Eq. ([7|, determined within this work for different LAB based scintillators. 
The upper limits for C are given for a confidence level of 95%. 



LAB admixture 


kB 


c 






[cm MeV-i] 


[cm^ MeV" 




2g/lPP0, 15mg/lbis-MSB 


0.0097 ±0.0002 


< 5.0x10" 


7 


2g/lPPO 


0.0098 ±0.0003 


<4.0 xlO" 


7 


3g/lPP0, 15mg/l bis-MSB 


0.0098 ±0.0003 


< 1.0 X 10" 


7 


3g/lPPO 


0.0094 ±0.0002 


<6.5 xlO" 


7 



Table IV. Detector properties assuming LAB with 2 g/1 PPO as target in SNO± and LAB with 3 g/1 PPO and 15 mg/1 bis-MSB 
in LENA. The SNO± target mass and energy resolution is the same as employed in [T^] for better comparison of the results 
and in agreement with ^Sj j. The LENA target mass and resolution is taken from [sj], where the expected resolution is quoted 
to be > 200 p. e. /MeV. The number of free protons Np is determined using a proton density in LAB of 6.31 x 10^^ cm""^ [30|. 
C is assumed to be zero. 



Detector Mass 


Np 


E resolution 


kB 


C 


[ktons] 


[10^^] 


{E in MeV) 


[cm/MeV] 


[cm^/MeV^; 


SNO± 0.8 


5.9 


5.0%/Vs 


0.0098 


0.0 


LENA 44 


325 


7.0%/Ve 


0.0098 


0.0 




Figure 1. Detected events sorted into a pulse-shape (PS) vs. pulse-height (PH) matrix 
in low gain (LG) mode and with selection of the energy of the incident neutrons. The 
structures labeled with 'n' and '7' are caused by charged particles resulting from neutron 
interactions with the scintillator and by Compton electrons from photon interactions. 
The photons result partly from ambient 7-radiation and partly from inelastic neutron 
scattering, n/7-separation, indicated by the solid line, can not fully separate neutron- 
and photon-induced events due to the small number of photoelectrons produced per 
scintillation light photon. The scintillator contained LAB ± 2 g/1 PPO + 15 mg/1 bis- 
MSB. 
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Figure 2. Pulse-height spectrum of Cs given in pulse-height channels. The calculated 
Compton electron spectrum before (blue solid line) and after (yellow shaded area) folding 
with the Gaussian resolution function is compared to the experimental spectrum (black 
filled circles) for LAB with 2g/l PPO and 15mg/l bis-MSB, measured in HG mode. 
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Figure 3. Relative pulse-height resolution AL/L for electrons as function of the light 
output L in electron-equivalent energy. The respective solute admixture to LAB is given 
in the legend. Shown are the data points acquired from source measurements in HG mode 
and the respective adapted resolution function Eq. ((3|. In the shown total uncertainties, 
the single contributions are added quadratically. 
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Figure 4. Position of the Compton edges in pulse-height channels as function of the 
maximum kinetic energy Ee of the Compton electrons. The data points are taken in 
HG mode. The linear fits describe the response of different LAB based scintillators to 
electrons. The respective solute admixture is given in the legend. In the shown total 
uncertainties, the single contributions are added quadratically. 




Figure 5. Calculated (red line) and measured (black histogram) proton recoil spectra at 
different incoming neutron energies. The measured spectra are extracted from the data 
sets for LAB with 2g/l PPO and 15mg/l bis-MSB. For the spectrum at 3MeV neutron 
energy, HG data is used, for the rest LG data. The calculation is performed using the 
NRESP7 code and does not consider a threshold (leading to counts starting at channel 1, 
where there is no data available) and photon induced background (leading to less counts 
at low pulse-heights above threshold compared to data). 
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Figure 6. Relative pulse-hight resolution AL/L for protons as function of the light output 
L in electron-equivalent energy for four different LAB based scintillators. The respective 
solute admixtures are given in the legend. In the shown total uncertainties, the single 
contributions are added quadratically. 
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Figure 7. Light output L, given in electron-equivalent energy, as function of kinetic 
energy E. Shown are the light output functions inserted in nresp7 for the initial (orange 
dashed dotted line) and the final (violet dotted line) proton spectrum calculations for 
LAB with 2 g/1 PPO and 15 mg/1 bis-MSB. 




Figure 8. Light output L, given in electron-equivalent energy, as function of kinetic 
energy E. Electron data (triangles) is fitted with Eq. ([5]) (grey dashed line) and proton 
data (circles) with Eq. ([7| (red solid line). The data is taken with LAB, 2g/l PPO and 
15 mg/1 bis-MSB. In the shown total uncertainties, the single contributions (listed in 
Tab. HI)) are added quadratically. For some values of L{E), the resulting error bars are 
smaller than the marker size. 




Figure 9. Same as Fig. [8] but for LAB with 2g/l PPO, without bis-MSB. 




Figure 10. Same as Fig.|H]but for LAB with 3g/l PPO and 15mg/l bis-MSB. 




Figure 11. Same as Fig.[8]but for LAB with 3g/l PPO, without bis-MSB. 




Figure 12. Proton quenching factor Qp as function of kinetic energy E for LAB based 
scintillators. The respective solute admixture is given in the legend. 
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Figure 13. Neutrino fluence on Earth induced by a lOkpc distant supernova, emitting 
a total energy of 3 x 10^"^ erg. The assumed mean neutrino energies are 12 MeV for !/e, 
15 MeV for Ve and 18 MeV for v^. 




1 2 3 4 5 

Proton Energy (MeV) 

Figure 14. True energy spectrum of recoiled protons in a ~0.8kton LAB detector (like 
SNO+) from f-p elastic scattering events, assuming the supernova neutrino fluences in 

Fig.na 
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Figure 15. Observed energy spectrum of recoiled protons in SNO+ considering proton 
quenching with kB — 0.0098 cmMeV"^ and an energy resolution of 5.0%/^/E. The su- 
pernova neutrino fluences in Fig. [TJ] are assumed. 
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Figure 16. Observed energy spectrum of recoiled protons in LENA considering proton 
quenching with kB — 0.0098 cmMeV"^ and an energy resolution of 7.0%/\^. The su- 
pernova neutrino fluences in Fig. [T3] are assumed. 




Figure 17. Observed energy spectrum of recoiled protons in SNO+ considering an energy 
resolution of 5.0%/\/!E and a supernova neutrino fluence as in Fig. 1131 Shown is the 
sum event spectrum for three different proton quenching scenarios, i.e. no quenching, 
kB = 0.0073 cmMeV"^ and kB = 0.0098 cmMeV"^ The red line indicates a threshold 
of 200 keV. 



